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Abstract. In this work we study the noise induced effects on the dynamics of short 
polymers crossing a potential barrier, in the presence of a metastable state. An 
improved version of the Rouse model for a flexible polymer has been adopted to mimic 
the molecular dynamics by taking into account both the interactions between adjacent 
monomers and introducing a Lennard-Jones potential between all beads. A bending 
recoil torque has also been included in our model. The polymer dynamics is simulated 
in a two-dimensional domain by numerically solving the Langevin equations of motion 
with a Gaussian uncorrelated noise. We find a nonmonotonic behaviour of the mean 
first passage time and the most probable translocation time, of the polymer centre 
of inertia, as a function of the polymer length at low noise intensity. We show how 
thermal fluctuations influence the motion of short polymers, by inducing two different 
regimes of translocation in the molecule transport dynamics. In this context, the role 
played by the length of the molecule in the translocation time is investigated. 



PACS numbers: 64.70.km, 83.10.Mj, 87.15.K- 



I Email: afiasconaro@gip.dft.unipa.it 
§ Email: spagnolo@unipa.it 




Noise driven translocation of short polymers in crowded solutions 



2 



1. Introduction 

The knowledge of the translocation dynamics of a molecule moving across a membrane 
or surmounting a barrier is a fundamental step towards the full comprehensions of 
the functioning of many biological systems. In the cell environment, DNA and 
RNA translocate across nuclear pores and many proteins work on the bases of their 
ability to go beyond a potential barrier. Reciprocal translocation between portions of 
chromosomes causes the genetic mutations which characterize several tumors [U E] . The 
cancer target therapy is based on a drug delivery mechanism that crucially depends on 
the translocation time of the chemotherapeutic molecules [SiSlE]. Polymer translocation 
has been found to play a key role in the acquisition of multi-drug resistance to 
cancer therapy [6]. The study of the transport of macromolecules across a nanopore 
is also important for technological applications. For example, the development of 
chemomechanical actuators [7J, the separation of molecules by liquid chromatography 
[H [9] or fast DNA sequencing techniques are both performed by forcing the molecule 
to move inside nano-channel devices. 

Many experimental studies on polymer translocation have been carried out in 
the wake of the pioneering work of Kasianowicz and collaborators (1996). In these 
experiments, single-stranded DNA (ssDNA) molecules are forced by a voltage bias to 
pass through an a-hemolysin (a-HL) pore imbedded into a membrane. The passage 
of the molecule inside the protein channel causes a reduction of the electrolyte ion 
current and the duration time of any current blockade is recorded. In this way, a 
linear relationship of the most probable crossing time Tp with the molecule length and 
an inverse proportionality law between Tp and the applied voltage were established 
[11] . The a-HL channel device has been extensively used for probing the transport 
dynamics of several types of polynucleotide (DNA or RNA) molecules having same 
length but different adenine and cytosine contents [iHl E]. A different nucleotide 
composition of the DNA brings about a significant change in the distributions of the 
duration times of the corresponding current drops. These findings confirm that the 
dynamics of biopolymer translocation across an a-HL channel is governed by pore- 
molecule interactions, which are dependent on the details of DNA sequences [15], the 
orientation and the driving voltage [IE]. Moreover, Tp scales as the inverse square of 
the temperature and the differences on the crossing time due to different nucleotide 
compositions are progressively reduced when the temperature increases [T3] . 

Recent advances in semiconductor technology have enabled the construction of 
synthetic nanopores for DNA translocation experiments [TTlfTS]. By using a silicon-oxide 
nanopore, it was found a power-law relationship between Tp and the polymer length, 
with an exponent of 1.27 [191 ISO]- However, a retarding effect on the crossing time has 
been observed as an effect of the voltage dragging itself, which could be stronger enough 
to disturb the biopolymer on finding the best orientation for the translocation [21]. In 
fact, other experiments show that DNA molecule's translocation speed in a solid-state 
nanopore is strongly dependent on electrolyte temperature, salt concentration, viscosity 
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and applied voltage bias [23j . 

In contrast to the proportionality law between the translocation time and the 
polymer length, experiments on the transport dynamics of DNA molecules driven inside 
an entropic trap array have shown longer crossing times for shorter molecules, suggesting 
the existence of a quasi-equilibrium state of the polymer during the passage [Ml EH] ■ 

The complex scenario of the translocation dynamics coming from experiments has 
been enriched by several theoretical and simulation studies |26l [271 l2Hl [291 HH [301 [211 
[321 [15]. A power law relationship between the translocation time of a long A^-segments 
chain molecule and the polymer length has been found, with an exponent of 2 in the 
presence of a free energy bias between the two sides of the membrane or 3 in the lack 
of adsorption [33]. The relationship between the crossing time and N becomes linear in 
the presence of strong interactions between the polymer and the pore walls p8] or in 
the case of an hairpin crossing mechanism of translocation pi] . 

In spite of the above contributions, the complicated biological environments and 
boundary conditions make the problem of polymer translocation still far from a common 
understanding and, in this framework, a detailed description of the transport dynamics 
of short chain molecules is missing. The motion of a polymer passing through a 
pore takes place in solutions where thermal fluctuations always affect the translocation 
dynamics. In this paper we present the results of our studies on the noise-induced 
effects on the transport dynamics of short polymers surmounting a potential barrier, in 
the presence of a metastable state. Molecular dynamics simulations are performed by 
modeling the polymer as a flexible chain molecule with harmonic interactions between 
adjacent monomers and a Lennard- Jones (LJ) potential between all beads. A bending 
recoil torque has also been included in our model. The polymer dynamics is simulated in 
a two-dimensional domain by numerically solving the Langevin equations of motion with 
a Gaussian uncorrelated noise. The dependence of the mean first passage time (MFPT) 
of the polymer centre of inertia from the length L of the chain molecule critically changes 
with the noise intensity [35]. In this paper we focus our research on short polymers, 
which usually show two different regimes of translocation depending on the ratio of L 
over the length of the channel [36] . We find a non-monotonic behaviour of both MFPT 
and function of polymer length at low noise intensity. In this context, the role 

played by the stiffness of the molecule in the translocation dynamics is also investigated. 
In section [2] we present our polymer chain model and give the details of the molecular 
dynamics simulations. The final results are described in section [3] and discussed in 
section [H 

2. Model and method 

2.1. The polymer chain model 

In our simulations, the polymer is modeled by a semi-fiexible chain of A^ beads connected 
by harmonic springs [37]. The contour length is defined as Lq = Nd, where d is the 
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equilibrium distance between adjacent monomers. Both excluded volume effect and 
van der Waals interactions between all beads are kept into account by introducing a 
Lennard- Jones potential. In this work, we investigate the dynamics of a linear chain 
molecule. In order to confer a suitable stiffness to the chain, a bending recoil torque is 
included in the model, with a rest angle = between two consecutive bonds. The 
polymer motion in the liquid solvent induces a velocity field which is felt by all the 
beads. To first order, we neglect this hydrodynamic effect. Such approximation implies 
that our model cannot be used to mimic the molecule behaviour in dilute polymeric 
solutions, but it appears to be much more appropriate for polymeric melts [38]. The 
potential energy of the modeled chain molecule is 

U = t/Har + f/Bcnd + t^LJ (l) 

where t/nar represents the energy required to extend the bond between two consecutive 
beads, [/Rend the work required to bend the chain, and [/lj the Lennard- Jones potential. 
Respectively, we have 

[/Har= 5^i^r(rM+l-rf)' (2) 

1=1 
N-1 

%end = ^ -^"6l(^j-l,i+l — ^o)^ (3) 

(4) 



i=2 



12 

a \ a 



t/Lj = 4eLj 

where K,. is the elastic constant, the distance between particles i and j, Kg the 
bending modulus, eLj the LJ energy depth and a the monomer diameter. 

2.2. Molecular dynamics simulations 

The effect of temperature fluctuations on the dynamics of a chain polymer escaping 
from a metastable state is studied in a two-dimensional domain. The polymer motion is 
modeled as a stochastic process of diffusion in the presence of a potential barrier having 
the form: 

UExt{x) = ax^ — bx^ (5) 

with parameters a = 3 ■ 10^'^ and b = 2 ■ 10^^. A three-dimensional view of f/gxt is 
plotted in figure [TJ The drift of the i^^ monomer of the chain molecule is described by 
the following overdamped Langevin equations: 

where U is the interaction potential, defined by equation ([1]), ^x{t) and iy{t) are 
white Gaussian noise modeling the temperature fluctuations, with the usual statistical 
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properties, namely {Ck(t)) = and {C,k{t)C,i(t + r)) = D5(k i)5{T) for (fc, / = x^y) and 7 is 
the friction coefficient. In this work we set the 7 parameter equal to one. The standard 
Lennard- Jones time scale is tlj = (mcr^/eLj)^^^, where m is the mass of the monomer. 
A bead of a single-stranded DNA is formed approximately by three nucleotide bases 
and then a ~ 1.5 nm and m k, 936 amu [15]. Orders of magnitude of the quantities 
here involved are nanometers for the characteristic lengths of the system (polymer and 
barrier extension) and microseconds for the time domain. Moreover, our simulation 
time ts is scaled with the friction parameter as ts = t/7, therefore we use arbitrary units 
for all computed translocation times. We set both and Kg equal to 10, eLj equal to 
0.1 and a equal to 3, in arbitrary units. The inter-beads rest length d is chosen equal 
to 5, while the number of monomers ranges from 4 to 40 units. 
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Figure 1. 3D-view and projection on z-x plane of the potential energy UExt, 
which is included in our system to simulate the presence of a barrier to be 
surmounted by the polymer. A sketch of the translocating chain molecule is 
shown. 

A sequence of 10^ numerical simulations has been performed for each of the 7 
different values of the noise intensity D and 11 different polymer lengths. The starting 
condition for all performed molecular dynamics simulations provides that the initial 
distance between two adjacent monomers of the chain is equal to the rest length d of 
the ideal spring connecting them. We have carried out our study by selecting an initial 
spatial distribution of the polymer with all monomers at the same x coordinate equal 
to Xo=0, which corresponds to the local minimum (metastable state) of the potential 
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energy of the barrier. Every simulation stops when the x coordinate of the center of 
mass of the chain reaches the final position at Xf=15. This value has been chosen after 
performing a set of simulations by assuming x/=30 and the noise intensities D = 0.3, 
1.0, 4.0 and 10.0, respectively. We have found that the polymer center of mass takes a 
very short time to travel from x=15 to x/=30 in comparison with the average escape 
time from Xi=0 to Xf=15. Therefore, by assuming Xf=15, we are confident that the 
probability for the polymer to be trapped back into the potential well is extremely 
low. The distribution of the translocation times is analyzed and the MFPT, Tp and the 
median of the distribution (hereafter 'the median') are calculated. 

3. Results 

The polymer center of mass crosses the line at x/ with translocation times having a 
distribution that crucially depends on the noise intensity D. The mean crossing time, 
the median and Tp are analyzed as a function of D for different values of the polymer 
length. In particular, figureEK shows a decreasing trend of the MFPT with increasing the 
intensity of the fluctuations affecting the monomers motion. This result is expected from 
the standard theory of the transport dynamics of a particle escaping from a metastable 
state [39]. However, when the noise intensity goes down D ~ 1, the polymer dynamics 
appears to be very sensitive to the molecule length, with shorter chains {N < 12) 
translocating more slowly (higher MFPTs) compared to longer molecules. In particular, 
the dynamics of polymers with iV ^ 10 and D ^ 1 shows interesting details inside the 
common behaviour, while the dependence of the MFPT on D converges to an almost 
unique path for greater values of A^. The crossing dynamics of chains with = 4 (blue 
line in figure [2^) resembles that of a single particle (see figure 2 in Ref. [40] and figure 1 in 
Ref. [S]), by showing a marked power law relationship between the MFPT and the noise 
intensity (straight line in a log-log plot). The comparison of the polymer behaviour with 
that of a single particle started outside the potential minimum, as studied in Refs. [lO] 
and [H], is possible because the exact initial point becomes irrelevant at strong noise 
intensities, where the power law behaviour is observed. 

Many experimental works report their results in terms of the most probable 
translocation time instead of the MFPT. For this reason, we have also investigated 
how the median and Tp depend on the noise intensity. The dependence of the median on 
D is very similar to that of MFPT, with a little exception aX D = 10. On the contrary, 
the diagram of Tp vs. D in figure [2]3 shows some diversities. The single-particle-like 
behaviour is still observed for = 4. However, when the polymer length increases, the 
index of the power law between Tp and D (the slope of the curve) progressively drops. 
For ^ 12 the translocation dynamics changes with respect to shorter chains and 
the dependence of Tp from the noise intensity becomes almost equal to that observed 
for the MFPT. In figure Et the ratio Tp/MFPT is plotted as a function of D. We 
can see that the two characteristic crossing times are not related to each other by a 
simple proportionality law for different polymer lengths and noise intensities. However, 
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Figure 2. Log-log plots of: a) Mean first passage time vs. noise intensity for 
11 different values of the polymer length; b) most probable translocation time 
vs. noise intensity; c) ratio of the most probable translocation time on MFPT 
vs. noise intensity. The colors correspond to different numbers N of chain 
beads: blue (iV = 4), red (iV = 6), green (iV = 8), yellow {N = 10), purple 
(A^ = 12), cyan {N = 15), grey (A^ = 20), orange {N = 25), black (A^ = 30). 



a common nonmonotonic behaviour of this ratio as a function of the noise intensity is 
observed. 

The different dynamics of translocation of short and long polymers is investigated 
as a function of the noise intensity. Figure [3] shows the dependence of the MFPT, the 
median and Tp on the polymer length (the number of chain beads) for four different 
values of noise intensity. In the diagram with D = 0.3 (figure [3^) two different regimes 
of translocation are clearly present: for A^ < 15 all characteristic times show a maximum 
at A^ ^ 7, while for A^ > 15 a common plateau is observed. This maximum is more 
evident in the MFPT than in Tp, but the height of both peaks is rapidly reduced when 
the noise intensity increases. This maximum occurs at a polymer length corresponding 
to the persistence length of the chain Lp, which represents the length over which the 
correlations in the direction of the tangent to the polymer are lost. For our discrete 
chain model, we have calculated the directional correlation: 

(cos(^(.,))) = ^^-T\^^-^)i^^^^-^) (8) 

as a function of the contour length sj corresponding to the j-th segment of the polymer 
chain, where Xk represents cos(6'(sfc)), 9{sk) is the angle between the tangent vector at 
one end of the chain molecule and the tangent vector at the contour length corresponding 
to the k-th segment of the polymer chain, x the mean of x. We have computed an 
ensemble average over 5 x 10^ different steric configurations of the polymer fluctuating 
in a fiat potential domain at the noise intensity D = 0.3 (the value for which the 
maximum of translocation time is observed in figure [3^). The directional correlation 
presents damped oscillations as a function of the polymer contour length, following 
an overall exponentially-decreasing trend. By best fitting this calculated correlation 
function with the expression (cos(6'(s))) = exp{—s/Lp), which is vahd for the continuous 
worm-like chain model, we have obtained Lp = 29.8. This value of Lp corresponds to 
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Figure 3. Mean first passage time (blue line), median of the crossing time 
distribution (red line) and most probable translocation time (green line) 
vs. polymer length (number of beads) for four different noise intensity D = 
0.3,1.0,4.0,10.0. The inset on each diagram shows an enlarged view on the 
y-scale of the region with > 20. 



AT ~ 6. 

In the short-length domain, the most probable translocation time shows an inversion 
of the dependence from the polymer length at higher values of D with respect to that 
observed at lower noise intensities. For D > 2 the most probable translocation time 
increases monotonically with the chain length and, for > 15, it reaches almost the 
same level of the MFPT. The inset of each diagram in figure [3] shows the fiat region 
of the corresponding plot in the range 20 ^ ^ 40 in an enlarged y-scaling. For any 
values of D, we observe a slightly increasing trend of the characteristic crossing times 
for longer molecules. 
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4. Conclusions 

In this paper we focus our study on the noise influence on the translocation dynamics 
of short polymers. We model the molecule as a chain of spatially-extended interacting 
beads. The transport dynamics is simulated in a noisy environment and in the presence 
of a potential barrier by solving the Langevin equation of motion for every single 
monomer. 

We flnd a maximum of the MFPT and 7~p dbS 8b function of the number of beads of the 
chain at low noise intensity. This maximum occurs at a polymer length corresponding to 
the persistence length of the chain Lp. This characteristic length represents a transition 
point between two different dynamical behaviours of the polymer approaching the 
barrier. The chain molecule having length lower than Lp moves like a flexible elastic rod 
and the addition of a monomer to the chain causes longer crossing times. For lengths 
greater than Lp the polymer can fold, assuming a snake-like shape. This conflguration 
helps the molecule into the translocation process, reducing the crossing times. For short 
molecules in little noisy environment, the relationship between Tp and the chain length 
is similar to that of MFPT. This means that both Tp and MFPT can be used without 
distinction to investigate the polymer translocation dynamics. On the contrary, when 
the thermal fluctuations are stronger, the MFPT and Tp show two distinct dependencies 
on A^. 

In order to compare our results with available experimental flndings, we restrict 
the following discussion on the behaviour of the most probable time. In a low noise 
environment, we conflrm that short polymers can travel more slowly than longer ones, 
as flrstly observed by Han and collaborators (1999). For greater values of polymer length 
(A^ > 25), our flndings are in agreement with the almost-linear relationship between the 
most probable translocation time and the polymer length [11] . At higher temperatures, 
the fluctuations dominate the transport dynamics of short polymers, canceling any 
dependence of the translocation time on the stiffness of the molecule. In this regime, 
longer molecules take longer times to cross the barrier. 

For polymer chains having length in a range close to Lp the noise acts as a trigger of 
the translocation dynamics. Low intensities of thermal fluctuations bring the molecule 
to cross the barrier by a transport mechanisms that is completely different from that 
observed at higher temperatures. The measured crossing times are a direct consequence 
of the noise influence on the molecule dynamical regime of translocation. Further work 
is required to include in our description a polymer-pore interaction term and to explore 
the noise induced modiflcations of the polymer dynamics in a sub-persistence length 
conflnement . 

The saturation of all the times represents an effect still not completely understood. 
A possible explanation could be related to the intrinsic stochastic behavior of the 
polymer dynamics: large number of monomers gives an average contribution to the 
motion of the polymer center of mass. When the number A^ increases, the fluctuations 
of the single monomer affect less and less the overall motion of the polymer. This 
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could be considered as a "coarse-grained" one-particle description in the limit — > oo. 
Moreover, the three characteristic translocation times saturates to almost the same 
value, indicating that the escaping time probability distributions tend to become more 
symmetric by increasing the number of monomers. A detailed investigation of these 
features will be the subject of a future work. 
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